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Abstract 


The  purpose  of  this  work  is  to  develop  a  large  signal 
signal  lumped  circuit  model  of  the  GaAs  MESFET  to  aid  in  the  designs 
of  microwave  MESFET  circuits.  The  circuit  elements  of  this  model 
are  obtained  either  directly  or  indirectly  from  the  DC  and  RF 
measurements  of  the  device  to  be  modeled. 

To  analyze  this  circuit  model,  a  nonlinear  circuit 
simulation  computer  program  is  written.  This  routine  is  base  on  a 
hybrid  time-frequency  domain  Iterative  algorithm  called  ‘multiple 
reflection  technique'.  To  Improve  the  speed  of  analysis,  an  accelerate 
convergence  scheme  is  Incoporated  into  the  multiple  reflection 
technique  for  the  first  time  to  analyze  three  terminal  device. 

The  validity  of  the  analysis  algorithm  is  first  checked  by 
comparing  the  simulation  results  of  a  MESFET  with  published  data. 
The  large  signal  model  developed  is  then  confirmed  by  comparing  the 
simulation  results  of  a  MESFET  modeled  in  this  work  to  the 
experimental  results. 
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Chapter  1  Introduction 


The  GaAs  metal  semiconductor  field  effect  transistor 
has  become  the  workhorse  device  in  the  microwave  industry.  The 
device  is  used  extensively  in  many  small-signal  applications  such  as 
low  noise  amplifier.  Due  to  the  exceptional  power  performance  of  the 
GaAs  MESFET,  which  is  comparable  or  superior  to  most  other 
commercial  microwave  solid-state  devices,  a  wide  range  of  power 
applications  have  been  seen.  The  applications  of  the  GaAs  MESFET 
such  as  for  power  amplifiers,  oscillators  and  mixers  are  growing 
rapidly. 

Despite  the  large  demands  for  GaAs  MESFET  as  power 
devices,  the  large-signal  design  techniques  for  these  circuits  are 
still  relatively  primitive.  This  situation  can  be  explained  by  two 
reasons:  1.  the  immaturity  of  the  GaAs  technology  in  general,  2.  The 
difficulty  in  characterizing  the  device  in  microwave  and  millimeter 
wave  frequencies.  Conventionally,  GaAs  MESFET  power  amplifiers  and 
oscillators  are  designed  by  using  small-signal  S-parameters  with 
tweaking  and  tuning  of  the  circuit  coming  later.  In  power  amplifier 
designs,  the  load-pull  method  can  also  be  used  to  determine  the 
optimum  load  impedance.  However,  this  load-pull  procedure  is 
laborious  and  cannot  predict  the  gain  performance  of  the  circuit. 

A  large  signal  model  of  the  GaAs  MESFET  can  aid  in  the 
design  of  a  nonlinear  MESFET  circuit  tremendously.  Several  research 
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efforts  have  tried  to  develop  a  numerical  model  of  the  GaAs 
MESFET[1][2][3],  These  numerical  models,  use  numerical  techniques 
such  as  a  finite  element  method  to  solve  the  nonlinear  differential 
equations  that  govern  the  device's  physical  properties.  These 
numerical  models  are  helpful  in  understanding  device  operation  and 
can  be  used  to  aid  in  the  design  of  the  MESFET.  However,  they  are 
impractical  for  circuit  level  design  purposes  due  to  the  long 
computation  time  required  for  these  models.  Mad  jar  and 
Rosenbaum[4],  Shur[5],  Chua  and  Sing[6]  developed  approximate 
analytical  models  for  the  GaAs  MESFET.  With  the  exception  of  Madjar 
and  Rosenbaum's  model,  none  of  these  analytical  models  have  shown 
validity  in  large  signal  circuit  design.  Furthermore,  this  type  of 
model  requires  the  detailed  knowledge  of  the  device's  physical 
parameters  which  often  cannot  be  obtained  accurately.  Willing, 
Raucher,  and  de  5ant1s[7]  developed  a  circuit  model  of  the  GaAs 
MESFET  by  measuring  the  device  S-parameters  at  various  bias  points. 
The  data  collected  is  then  used  to  extract  linear  element  values  of 
the  circuit  model.  The  voltage  dependencies  of  the  nonlinear 
elements  are  obtained  by  using  a  polynomial  curve-fitting  to  the 
S-parameters  measured  at  each  bias  point.  Although  this  approach 
was  shown  accurate,  it  has  the  drawback  of  being  very  tedious. 
Tajima,  et  a1.[8][9],  Materka  and  Kacprzak[10)  used  a  quasi-static 
approach  to  develop  their  circuit  model.  In  this  approach,  the  voltage 
dependencies  of  nonlinear  elements  are  assumed  to  be  the  same  as  in 
the  DC  condition.  They  used  DC  current-voltage  characteristics  of  the 
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ME5FET  to  establish  the  voltage  dependencies  of  the  nonlinear 
elements  and  used  S-parameter  optimization  to  obtain  linear  element 
values.  Their  results  showed  this  approach  is  accurate  at  least  up  to 
18  GHZ.  In  this  study,  a  crcuit  model  similar  to  the  one  used  by 
Materka  and  Kacprzak  is  employed.  A  systematic,  yet  simple 
technique  of  identifying  circuit  element  and  circuit  parameter  values 
is  developed. 

The  circuit  model  thus  developed  needs  to  be  analyzed  by 
a  nonlinear  circuit  analysis  technique.  Frequently  used  nonlinear 
circuit  analysis  methods  can  be  divided  into  two  classes:  i.  the  time 
domain  analysis,  2.  the  hybrid  analysis  that  iterates  between 
frequency  domain  and  time  domain. 

The  time  domain  approachs  includes  the  direct 
integration  method  and  the  shooting  method.  The  direct  integration 
method  is  a  brute  force  method  which  solves  the  nonlinear 
differential  equations  that  describes  the  nonlinear  circuit  at  each 
discretized  time  step  until  it  reaches  steady  state.  This  approach 
often  requires  long  computation  time  if  the  transient  time  is  long.  To 
avoid  long  transient  analysis,  Colon  and  Trick[l  1]  proposed  a  shooting 
method  that  uses  Newton's  iteration  algorithm  to  find  a  steady  state 
solution.  Skelboet  1 2]  used  a  different  kind  of  shooting  method  with 
an  extrapolation  algorithm.  Even  with  the  shooting  method,  The 
computation  time  for  the  time  domain  approach  is  long  compared  with 
a  hybrid  type  of  analysis. 

Nakhla  and  Vlach  used  a  hybrid  analysis  called  ’harmonic 


balance  method’.  In  this  algorithm,  the  circuit  under  analysis  is 
decomposed  into  two  parts,  the  nonlinear  subnetwork  and  the  linear 
subnetwork.  The  linear  subnetwork  is  described  by  linear  equations 
in  the  frequency  domain  and  the  nonlinear  network  is  described  by 
nonlinear  differential  equations  in  the  time  domain.  Interconnecting 
the  linear  and  nonlinear  network  can  be  achieved  by  relating  the  time 
and  frequency  domain  by  the  Fourier  transformation.  The  harmonic 
compoments  of  voltages  an/or  currents  at  the  interconnection  of  the 
linear  and  nonlinear  subnetworks  are  chosen  as  an  independent 
variable  of  the  system.  These  variables  are  optimized  by  an  iterative 
algorithm,  so  that  the  linear  equations  of  the  linear  network  and  the 
nonlinear  differential  equation  of  the  nonlinear  network  are  all 
satisfied.  With  a  suitable  optimizing  method,  the  harmonic  balance 
method  can  be  quite  efficient.  However,  the  success  of  this  approach 
often  depends  on  good  initial  guess  values  of  the  voltage  and/or 
current  harmonic  components.  If  the  initial  guess  values  are  not 
close  enough  to  the  actual  solution,  the  solution  might  converge  to 
local  minimums  or  not  converge  at  all. 

Kerri  14],  proposed  a  different  kind  of  hybrid  technique 
which  is  called  the  'multiple  reflection  method’  This  method  also 
iterates  between  the  time  and  frequency  domain.  It  requires  no 
initial  guess.  However,  the  speed  of  convergence  is  relatively  slow 
In  the  present  work,  the  multiple  reflection  technique  is  modified  by 
adding  a  ’voltage  update  algorithm'.  For  MESFET  circuit  analsyis,  the 
modified  multiple  reflection  method  is  found  to  be  much  more 
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efficient  than  the  original  multiple  reflection  method 

Chapter  2  of  this  thesis  describes  the  GaAs  MESFET 
model.  The  large  signal  model  is  evolved  from  the  conventional  small 
signal  lumped  circuit  model.  From  an  initial  sensitivity  study,  some 
of  the  circuit  elements  are  identified  as  nonlinear  elements.  By 
means  of  the  quasi-static  approach,  these  elements  are  modeled  by 
nonlinear  equations.  The  linear  element  values  are  obtained  from  DC 
and  RF  measurements. 

Chapter  3  is  devoted  to  the  large  signal  curcuit  analysis 
techniques.  The  multiple  reflection  method  is  introduced  first.  In 
this  method,  the  FET  circuit  is  divided  into  a  nonlinear  part  and  linear 
part.  The  solution  is  arrived  at  by  solving  the  nonlinear  and  linear 
circuit  iteratively.  A  voltage  update  procedure  is  then  added  to  the 
multiple  reflection  method  to  increase  convergence  speed. 

Chapter  A  presents  the  simulated  results  and  the 
experimental  data.  The  validity  of  the  modified  multiple  reflection 
method  is  first  confirmed  by  comparing  the  data  obtained  with 
Materka's  MESFET  circuit  simulations.  An  arbitrarily  chosen  device  is 
modeled  using  the  techniques  in  chapter  2.  Simulated  results  and 
measured  data  of  this  device  are  compared. 

Chapter  5  forms  the  conclusion  and  suggests  future 
research  paths.  Possible  related  future  research  topics  can  be  the 
modeling  of  a  dual  gate  MESFET,  the  analysis  of  other  large  signal 
devices  such  as  FET  mixer,  FET  oscillator  and  FET  multiplier. 


Chapter  2  GaAs  MESFET  Modeling 


The  purpose  of  this  research  is  to  develop  a  practical 
model  of  the  GaAs  MESFET  to  aid  in  the  large  signal  GaAs  MESFET 
circuit  lesign.  This  model  must  meet  four  requirements. 

1.  The  modeling  technique  should  be  general,  so  that  it  applies  to 
devices  with  various  dimensions  and  physical  parameters. 

2.  The  procedures  of  extracting  circuit  parameters  should  be 
reasonably  simple. 

3.  The  model  can  be  analyzed  in  a  computationally  efficient  manner. 
A.  The  model  should  be  reasonably  accurate  in  predicting  device 

performance. 

2.1  The  large  signal  equivalent  circuit 

The  large-signal  equivalent  circuit  is  derived  from  the 
conventional  small-signal  equivalent  circuit.  Figure  2.1  is  such  a 
small-signal  equivalent  circuit.  The  position  of  each  lumped  element 
in  an  actual  device  is  illustrated  in  figure  2.2.  Ld,  Lg  and  Ls  are  the 

bound  wire  inductances.  Rg  is  the  metallization  resistance  at  gate. 
Rd  and  Rs  consist  of  the  ohmic  contact  resistances  and  the 
resistances  in  the  active  layer.  Rj  is  the  charging  resistance  to  the 
gate-source  capacitance.  Cds  is  the  capacitance  that  exists  in  the 


channel  region.  R0  Is  the  output  resistance.  Precisely  speaking,  all 

these  elements  exhibit  certain  degrees  of  nonlinearity.  By  measuring 
the  S-parameters  of  the  device  at  sufficiently  different  bias  points, 
and  curve-fitting  the  small  signal  equvalent  circuit  to  the 
S-paramenters  measured  at  each  bias  point,  the  element  values  of  the 
small-signal  equvalent  circuit  versus  bias  voltages  curve  can  be 
established.  It  was  found  that  Ls,  Lg  Ld,  Rd,  Rs,  Rg  and  Rj  do  not 

exhibit  strong  bias  dependency[7).  These  elements  are  considered  as 
linear  elements  in  this  study.  At  this  point,  a  priliminary  large  signal 
equivalent  circuit  can  be  drawn(fig.  2.3).  Notice  that  R0  and  Gm  are 

combined  to  a  single  nonlinear  element  lch(Vg,  Vd).  Taking  into 

account  gate  voltage  swing  to  positive  voltage  and  the  drain 
break-down  phenomenon,  forward  gate  bias  current  and  drain  to 

gate  current  lbr  are  added  to  form  the  final  large  signal  equivalent 
circu1t(f1g.  2.4). 

12.  Extracting  linear  element  values 

Rg,  Rd,  Rs,  Lg,  Ld  and  Lg  are  identified  as  linear 

elements.  Their  values  are  assumed  to  be  voltage  independent.  In 
general,  these  elements  have  small  values  and  are  difficult  if  not 
impossible  to  measure  directly.  Often,  these  values  are  established 
by  measuring  device  5-parameters  at  certain  bias  conditions  and 
curve-fitting  the  small-signal  equvalent  circuit  to  the  5-parameters 


As  Diamond  and  Lavironl  1 5]  have  suggested,  the  S-parameter 
measurements  of  a  device  at  VDS-  0.  V,  permits  a  more  accurate 

evaluation  of  device  parasltics  because  the  device  can  be  represented 
by  a  simpler  equivalent  circuit.  Figure  2.5  is  the  equivalent  circuit  at 
Vqs  -0.  The  channel  region  under  the  gate  is  depleted  evenly  due  to 

the  lack  of  drain  bias  voltage.  This  depleted  channel  region  is 
modeled  by  the  distributed  R-C  network  in  the  equivalent  curcuit. 
Curtice  and  Camisa[16]  found  that  when  trying  to  optimize  the  circuit 
element  values  of  figure  2.5  to  curve  fit  the  measured  S-Parameters 
at  VDS-0,  with  small  difference  in  the  error  function,  the  optimum 

values  of  Rg,  Rg  and  Rd  vary  widely  depending  upon  the  optimization 
method  and  the  initial  guess  values.  This  suggests  the  parastic 
resistance  values  Rg,  Rd  and  R$  obtained  in  this  way  may  not  be 

reliable.  In  this  work,  a  DC  measurement  technique  called  Fukui's 
method!  17]  will  be  used  to  determine  parasitic  resistance  values. 
The  unknown  element  values  of  figure  2.5  are  now  reduced  to  6(R, 
C,CdS,  Lg,  Ld  and  Lg).  Now,  the  bound  wire  inductance  values  Lg,  Ld 

and  Ls  can  be  optimized  with  greater  accuracy. 

A  Hughes'  PF-6000  GaAs  MESFET  is  chosen  to  be  analyzed. 
The  device  has  a  gate  length  of  0.6jim  and  a  gate  width  of  960  jim. 
The  following  procedure  are  taken  to  evaluate  the  linear  element 
values. 

Since  the  parasitic  resistance  values  are  small  (on  the 


order  of  1  ohm),  the  measurement  system  needs  to  be  carefully 
calibrated.  The  Rg+Rs  value  can  be  approximated  by  measuring  the 
slope  of  VgS  vs.  lgS  at  high  current  level  with  the  drain  at  ground  and 

the  source  open.  Similarly,  Rg+Rg  can  be  calculatedffig.  2.6).  The 

Rg+Rg  value  is  extrapolated  by  the  following  steps: 

a.  measuring  the  drain  l-V  curve  in  the  linear  region  (figure 

2.7). 

b.  plotting  the  slope  of  V^g/lgg  at  each  gate  bias  voltage  vs. 
X,  where  X  is  defined  as, 

X={  1  -SQRTK  Vb-Vgs)/(  Vp+  Vb)ir 1  (2.1) 

Vg  is  the  built-in  voltage  for  the  Schottky-gate  and  Vp  is 
the  pinch-off  voltage. 

The  linear  extrapolation  of  this  curve  to  the  Y  axis  gives  the  value  of 
Rs*Rd  (figure  2.8).  From  the  values  of  Rg*Rd,  Rg+Rs  and  Rg+Rg,  the 

values  of  Rg,  Rg  and  Rs  can  be  subsquently  calculated. 

With  Rg,  Rg  and  Rg  known.  Lg.  Ld  and  Lg  can  be  optimized 

to  the  S-parameters  measured.  This  is  conveniently  done  by  using 
Super-Compact!  1 8]. 


2.3  Modeling  of  nonlinear  elements 

Nonlinear  elements  are  the  ones  that  exhibit  bias  voltage 
dependent  values.  The  quasi-static  approach  is  taken  to  model  the 
nonlinear  elements.  In  this  approach,  the  assumption  is  made  that  the 
voltage  dependency  of  the  nonlinear  element  at  high  frequency  can  be 
determined  from  its  dependency  on  DQ  bias  voltage 

a.  Channel  current  and  breakdown  current 

The  nonlinear  expression  for  channel  current  and  drain  to 
gate  breakdown  current  are  determined  by  the  DC  drain  current  versus 
the  drain  and  gate  bias  voltages  curves.  Figure  2.9  is  the  measured 
curve  of  the  Hughes  PF-6000  FET.  The  drain  current  is  equal  to  the 
sum  of  the  channel  current  and  the  breakdown  current. 


lds3,ch+,br 


(2.2) 


The  breakdown  current  lbr  is  significant  only  at  high  Vds  and  high  VgS 

values,  where  the  drain  to  gate  voltage  V^g  is  high.  The  negative 

differential  rsistance  observed  at  high  gate  voltage  after  saturation 
is  due  to  the  Gunn-domain  effect  which  is  caused  by  the  velocity 
saturation  characteristics  of  GaAs[l9).  The  Gunn-domain  effect 
produces  a  dipole  layer  in  the  active  layer  near  the  drain  side  of  the 
gate.  This  dipole  layer  also  produces  capacitance  effects  and  will  be 
discussed  later.  This  negative  resistance  phenomenon  is  especially 


pronounced  in  a  short  channel  GaAs  FET,  as  in  the  case  of  a  microwave 
GaAs  FET.  In  short  channel  devices,  a  high  electric  field  region  exists 
in  the  channel.  Electrons  that  enter  this  high  field  region  are 
accelerated  to  a  velocity  as  high  as  twice  the  equilibrium  velocity 
before  they  relax  to  the  equilibrium  veloci ty[20]  This  velocity 
overshoot  shifts  the  dipole  layer  into  a  gap  between  the  gate  and 
drain.  Figure  2.10  shows  the  location  of  this  dipole  layer 

the  emperical  expression  proposed  by  Naterka[10]  is 
modified  to  model  the  channel  current  lch, 

*ch=,dss  x  ^ ~Vg(t-f)/Vp)^ 

xtanh[cVd/(Vg(t-r)-Vp)]  Vp<Vsat 

=ldssx(|-Vg  (t-o/y2  (23> 

X  tanh[cvd/<vg(t-  r)-Vp)] 

-Vdxg0/(Vs-Vg(t-t))<I  vd>Vs3t 

where  V  =V  ♦  y  v 
p  po  d 

The  breakdown  current  is  modeled  by  an  exponential  expression, 
'br=lsr  exP(  P  Vpg)  (2.4) 

It  should  be  noted  that  even  though  the  expression  for  l^r  is  similar 
to  the  one  of  diode  forward  current,  they  represent  different  physical 
phenomenons.  The  parameters  ldss,  c,  g0  Vpo,  y,  Vs,  Vsat,  !sr  and  p 


are  optimized  to  the  measured  l-V  curves.  Using  these  two  emperical 
expression,  the  simulated  l-V  curves  for  the  PF-6000  MESFET  is 
shown  In  figure  2. 1 1 .  The  parameter  x  in  (2.3)  is  the  propagation  delay 
time  due  to  the  finite  velocity  of  the  charge  carriers  which  travel  in 
the  channel  region.  Electrons  entering  the  channel  at  the  source  end 
travel  for  a  time  x  before  it  reaches  the  drain  end.  This  propagation 
delay  time  is  a  basic  limitation  of  the  device.  Typically,  t  is  on  the 
order  of  lOps.  At  low  operating  frequencies,  xcan  be  neglected.  As 
frequencies  get  higher,  x  becomes  an  increasingly  important  factor 
To  calculate  t,  the  velocity  destribution  ll(x)  of  carriers  along  the 
channel  region  needs  to  be  known,  and 

o 

T«  Jl/U(x)dx  (2.5) 

s 

Since  D(x)  is  bias  dependent,  x  is  also  bias  dependent.  Due  to  the  high 
field  effect,  the  velocity  distribution  D(x)  cannot  be  easily 
calculated.  For  simplicity,  x  Is  assummed  constant  in  this  study. 
This  approximation  should  be  reasonable  in  the  saturation  region  of 
the  l-V  curves  where  the  charge  carriers  travel  at  the  saturation 
velocity  in  most  of  the  channel  region.  This  constant  value  of  xcan  be 
conveniently  obtained  by  using  the  small-signal  equivalent  circuit  of 
the  MESFET  in  the  Super-Compact,  which  has  a  time  delay  option.  In 
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Figure  2.12 


voltage  dependence  of  f.,,s 


this  equivalent  circuit,  the  value  of  t  is  optimized  to  fit  measured 
5-parameter  values  at  a  normaly  biased  condition  (eg.  VG5=-2.V, 

VD5=7V  for  the  PF-6000  MESFET). 

b.  Forward  gate-bias  current 

The  forward  gate  bias  current  is  modeled  by  the  diode 
forward  bias  current, 

lf=lso  exp(«.Vg)  (2.6) 

lso  and  «.  can  be  obtained  by  measuring  versus  VgS  with  the  drain 
terminal  open  ,and  plotting  vs.  Vg  on  log  paper,  where 

Vvgs-,f(Rs+Rg)  (27) 

The  slope  of  this  curve  equals  a,  and  the  intercept  on  the  Y-axis 
equals  lso. 

c.  Gate  to  source  and  gate  to  drain  capacitances 

The  gate  to  source  capacitance  CgS  consists  of  2 

components: 

cgs=csp+csd  (2.8) 

Csp  is  the  capacitance  due  to  gate  and  source  contact  metalization. 

CS(j  is  the  gate  to  source  depletion  capacitance.  Similarly,  Cgd-the 
gate  to  drain  capacitance  is  expressed  as 


Cdp  is  the  capactance  due  to  gate  and  drain  contact  metalization.  Cdd 
is  the  gate  to  drain  depletion  capacitance.  Cgd  and  Cdd  together 
contribute  to  the  total  gate  depletion  capacitance.  Willing,  et  al., 
have  shown  that  CgS  and  Cgd  are  both  gate  voltage  and  drain  voltage 

dependent[7],  with  CgS  being  weakly  drain  voltage  dependent.  To 

model  CgS,  the  depletion  capacitance  Csd  is  assumed  to  be  gate 

voltage  dependent  only  and  is  represented  by  the  Schottky  diode 
depletion  capacitance  expression, 

csd(  vg)=cso/SQRT(  1  -  V' V  for  V  8Vb  (2  1 0) 

where  Vb  Is  the  build-in  voltage  of  the  Schottky  barrier.  For  Vg>.8Vb, 

the  charge  storage  capacitance  is  assumed  to  dominate  and  increases 
linearly  with  Vg  with  Its  slope  equal  to  the  slope  at  Vg=.8Vb.  The 

metalization  capacitance  is  constant  and  can  be  calculated  by  using 
the  closed  form  formula  for  coupled  microstrip  lines[2l].  Now,  the 
total  expression  for  Cgs  can  be  represented  as 
CgS=CS0/5QRT(  1  -Vg/ VbKsp  ,for  Vg<.8Vb  (2 . , , ) 

-2.236CS0*  5.59Cgo/Vb  x  (Vg-.8Vb)+Csp  ,  for  Vg>.8Vb 
Figure  2. 12  shows  this  gate  voltage  dependence  of  CgS. 

The  gate  to  drain  cpapcitance  Cgd  decreases  rapidly  as 
Vfj  increases  in  the  linear  region  of  the  l-V  curves.  After  saturation, 


the  shape  of  the  depletion  region  under  the  gate  at  the  drain  side  does 
not  change  much  with  an  increase  in  drain  voltage  and  the  depletion 
capacitance  Cdd  becomes  small.  In  this  model,  Cdd  is  assumed  to 

decrease  linearly  to  0.  at  the  saturation  voltage  Vsat.  The  gate 
voltage  dependence  of  Cdd  is  the  same  as  that  of  Csd  in  this  region. 
The  total  gate  to  drain  capacitance  Cgd  is 
Cg<rCgo/SQRT(  I -Vg/V„)  X  ( 1  - Va/Vsat>.  Cdp  ,  for 

”cdp  *  ^or  vd>vsat 

(2.12) 

Cdp  the  metallization  capacitance  can  be  calculated  again  by  closed 

form  expression  for  coupled  mtcrostrlp  lines.  The  curve  for  Cgd  is 

shown  in  figure  2.13.  Cgo  of  (2.12)  and  Cgo  of  (2.1 1)  can  be  obtained 

by  drawing  the  device  equivalent  circuit  at  a  zero  biased  condition  as 
shown  in  figure  2.14,  and  by  using  the  5-parameters  measured  at  this 
bias  condition  to  optimize  Cgo  and  Cso  values. 

d.  Rj  and  Cd8 

Rj  is  the  gate  to  source  capacitance  charging  resistance. 
The  product  of  R|  and  Cgs  is  the  gate  to  source  capacitance  charging 
time  constant.  Since  Cgs  is  bias  dependent,  Rj  is  also  bias  dependent. 


/  f 


The  magnitude  o(  Rj  Is  usually  small.  For  simplicity,  Rj  Is  assumed 
constant  in  the  model.  The  value  of  Rj  can  be  optimized  by  using 
Super-Compact. 

Cds  is  the  capacitance  that  exits  in  the  channel  region 

Its  value  is  expected  to  be  small  before  saturation.  After  saturation, 
its  value  increases  due  to  the  formation  of  the  Gunn  domain  dipole 
layer.  For  the  PF-6000  MESFET,  Cds  is  found  to  be  .2pf  at  zero  biased 

condition(eg.  VD5=0.,  VG5*0.)  and  is  ,4pf  at  VG5=-2.V,  VD5=7.V.  In 

this  model  a  constant  Cds  value  is  assumed.  Willing,  et  al.[7j,  has 

stated  this  approximation  is  valid  at  least  in  the  saturation  region. 
The  value  of  Cds  can  be  optimized  at  a  bias  point  in  the  saturation 

region. 


Chapter  3  Large  Signal  Analysts 


3.1  Introduction 

A  nonlinear  circuit  analysis  computer  program  is  needed 
to  simulate  the  large  signal  model.  The  multiple  reflection  method  is 
chosen  to  analyze  the  nonlinear  circuit.  The  advantage  of  this  method 
over  the  popular  harmonic  balance  method  is  in  its  reliable 
convergence  nature  and  no  initial  guess  requirement.  The 
disadvantage  is  the  longer  computation  time  required.  At  the  initial 
development  stage  of  the  MESFET  model,  reliability  is  certainly  more 
important  than  efficiency  of  the  simulation  method.  In  this  study,  the 
multiple  reflection  method  is  improved  by  using  an  accelerated 
convergence  scheme.  This  modified  multiple  reflection  technique  is 
found  to  be  more  efficient  than  the  original  one  in  this  study. 

In  the  multiple  reflection  algorithm  and  other  hybrid 
time-frequency  domain  type  methods,  the  circuit  under  analysis  is 
divided  into  2  parts:  the  FET  intrinsic  circuit  and  the  external 
embedding  circuit.  The  external  embedding  circuit  consists  of  an 
input  matching  circuit  and  an  output  matching  circuit.  The  intrinsic 
FET  circuit  is  characterized  by  a  set  of  nonlinear  differential 
equations  in  the  time  domain.  The  external  embedding  circuit  is 
conveniently  represented  by  a  set  of  linear  phasor  equations  in  the 
frequency  domain.  These  2  sets  of  equations  must  be  solved 
simultaneously.  In  the  harmonic  balance  method,  the  harmonic 


compments  of  the  voltage  at  the  Junctions  of  the  FET  circuit  and 
external  circuit  are  chosen  as  independent  variables.  These  variables 
are  optimized  using  iterative  methods  such  as  the  steepest  descents 
method  and  the  Newton-Raiphson  method.  The  solution  is  arrived  at 
when  the  RMS  difference  of  the  currents  at  the  junction,  resulting 
from  each  set  of  equations  are  minimum.  With  a  fast  optimization 
routine,  this  harmonic  balance  approach  can  be  quite  efficient. 
However,  there  are  two  major  problems  with  this  approach.  First  of 
all,  it  requires  an  initial  guess.  The  speed  or  even  the  success  of  this 
approach  depends  strongly  on  how  good  the  initial  guess  is.  Second, 
the  optimization  routine  might  converges  to  local  minimums  which 
makes  the  harmonic  balance  mehtod  unreliable. 

3.2  The  Multiple  reflection  technique 

In  the  multiple  reflection  technique,  two  ficticious 
transmission  lines  are  assumed  to  exist  between  the  intrinsic  FET 
circuit  and  the  input-output  matching  networks,  as  shown  in  fugure 
3.1.  The  transmission  lines  are  assumed  lossless  and  are  large 
integer  multiples  of  wavelength  long.  With  this  assumption,  the 
steady  state  solution  is  not  changed  by  inserting  these  transmission 
lines. 

At  t-O,  the  gate  and  drain  input  voltages  are  applied  to 
the  transmission  lines.  The  initial  incident  voltage  wave  to  the  FET 
circuit  from  the  gate  side  is: 


(3.1) 


Vtj1(t)-Vgg  Z0/SQRT(Z02^Zg,2)  X  cos(Mt-phase) 

+VBBXZ0/(Zo*Zgo) 

.  where  phase*tan~ 1  [lm(Zg j  )/(Z0*Re(Zs , ))) 

The  initial  incident  voltage  from  the  drain  side  is 

vLi'(t>*VDDXZo/<VZLo>  (32> 

The  superscript  refers  to  the  iteration  number.  These  incident  waves 
reach  the  FET,  where  the  FET/ transmission  line  interface  can  be 
described  by  the  circuit  of  figure  3.2.  This  circuit  is  analyzed  in  time 
domain  and  will  be  discussed  in  next  section.  After  the  FET  circuit 
reaches  steady  state,  the  reflected  voltage  waves  are  calculated  as 

vlr,«)-vgsl<t,'lgl(t)Zo 

vLr'<t>-vds'«>-|d'«>zo  <33> 

These  voltage  waves  can  be  represented  by  Fourier  series, 

•o 

Vfrhthl  Vjr'(n©)  ej0®* 

n-0  (3.4) 


VLrht>J.  V^no^ei"®1 

n=0 

Let  the  operation  of  calculating  Vjr(m©)  and  VLr(m®)  be  F|t  J  and 
flU  respectively,  then 
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Vjr(na>)=F,(V1j(t)JVLi(t)l 

VLr(nto)=FL(V11(t),Vu(t)) 
n-0.1.2,. ... 


(3.5) 


These  reflected  waves  travel  back  to  the  embedding  circuit  and  are 
reflected  again. 

The  new  incident  waves  are  now  equal  to  these  reflected  waves  plus 
the  initial  incident  waves 


Vit2(t).V^(t)^  Vjpkn®) r.(n»)e-Jna,t 


(  3.6) 


VLi2(t)=Vu1(t)-2  VLr  1  (na>>r  (ntoJeJH0®1 


TjCnofl)  and  rL(n»)  are  the  reflection  coefficients  at  the  interfaces  of 
the  ficticious  transmission  lines  and  the  input/output  embedding 
circuits,  respectively,  r^nw)  andr^nw)  are  calculated  as 

rj(n®)”(Zsn-Z0)/(Zsn+Z0) 


rL(n «WZLn-Z0)/(ZLn*Z0> 
n-0, 1,2,3,... 


(3.7) 


*  '■-*  1  -g  \-v  7T  ^  'jv  1  «V  ^  ’■  V  ■  ’V  *»'  *>."  *  W  ^L"  v  ffvTv/  HT  v*  H."  K."  V.*” ' v^v^' v- j  r  vr  -.-■  v-/  v; 


This  process  continues  until  the  solution  converges  to  a  steady  state 
value.  The  multiple  reflection  scheme  is  illustrated  in  figure  3.3 
Convergence  occurs  when  the  RMS  difference  of  the  reflected  waves 
on  two  consecutive  iterations,  is  less  than  a  tolerance  value. 

The  value  of  the  characteristic  impedance  of  the 
ficticious  transmission  lines  will  affect  the  convergence  speed.  Z0 

should  be  chosen  close  to  the  value  of  ZS|  in  order  to  minimize  the 

reflection  coefficient  at  the  input  end.  Z0  also  affects  the  time  for 

the  intrinsic  FET  circuit  to  reach  steady-state  and  therefore  should 
not  be  too  large.  The  values  for  characteristic  impedance  of  the  input 
transmission  line  and  output  transmission  line  need  not  be  the  same. 

During  iterations,  the  gate  or  drain  voltage  may  exceed 
the  limiting  value  of  the  model.  This  can  be  prevented  by  setting  the 
DC  source  and  load  resistances  to  z0,  even  though  in  practice  they 

have  small  values.  This  change  has  no  effect  on  the  gate  circuit 
because  the  gate  current  has  a  very  small  DC  component.  The  DC 
voltage  drop  in  the  ficticious  resistance  can  be  compensated  in  the 
drain  by  offsetting  the  drain  bias  voltage. 

3.3  Time  domain  analysis  of  the  intrinsic  FET  circuit 

During  iterations,  the  circuit  of  figure  3.2  needs  to  be 
analyzed  in  the  time  domain.  The  inductance  Lg  and  lg  are  linear 

elements  and  can  be  moved  outside  the  intrinsic  FET  circuit.  Let 
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Rl“zo+Rg  and  R2aZo+Rd  The  c,rcu1t  of  figure  3.2  is  written  in  the 

form  of  figure  3  4  To  analyze  the  circuit  of  figure  3  4,  3  loops 
equations  can  be  written 

Vg(t)- vp(t)- V0(t)+R|  l  (t)-0  (3  8) 

-2Vu<t)*(R , ♦R2)lp(t)> Vp(t)*Rj  lg(t)-R2  ld(t)*2Vu(t)-0.  (3.9) 

2Vii(t)+R2lp(t)“Rslg(t)_(R2+Rs)ld(t)"vd(t)'Ls(ld(t)+l9(t))=0  (3  10) 
and  three  nonlinear  equations  can  be  written 

lg(t)=Cgs(Vg(t))  Vg(t)  ♦  lso  exp(ot  Vg(t))  (311) 

lp(t)=Cgd(Vg(t),  Vd(t))Vp(tHsrexp(-p  Vp(t))  (3.12) 

ld(t)=lch(Vg(t-r),Vd(t))+Cds  Vd(t)  (3.13) 

These  set  of  6  equations  can  be  solved  by  using  numerical  integration 
methods.  The  initial  condition  for  voltages  and  currents  of  the  first 
iteration  are  obtained  by  performing  a  DC  analysis  of  the  circuit  in 
other  words,  it  is  assumed  DC  voltages  have  been  applied  to  the 
circuit  for  t<0. 

3.4  Modified  multiple  reflection  technique 

Hicks  and  Khan  [23]  proposed  a  voltage  update  algorithm 
to  analyze  nonlinear  circuits.  In  this  method  the  circuit  under 
analysis  is  decomposed  to  two  parts(figure  3.5).  The  linear  portion 


and  the  nonlinear  portion.  The  algorithm  is  described  in  the  following 
steps 

1  Estimate  the  harmonic  components  of  the  periodic  voltage 
at  the  interface  V^tmtn)  (superscript  indicated  iteration 
number). 

60 

2.  Apply  Vn  1  (0=1  Vn  1  (rntatej^*  to  the  nonlinear  circuit,  and 

m=1 

calculate  ln  1  (t). 

3.  lL1(t)=-ln](t),  take  the  Fourier  transform  of  lL  1  (t)  to 
frequency  domain. 

4  Apply  ^(mco)  to  the  linear  network  and  calculate  ^'(mw) 
for  m=0,l,2 _ 

5.  New  harmonic  components  of  the  voltage  applied  to  the 
nonlinear  network  are  updated  as, 

Vnk+  1(mco)=Vnk(mci))+Pm  (VLk(m»)-Vnk(ma>))  ,0<Pm<1  (3  14) 

The  process  continues  until  the  RMS  difference  of  Vn(t)  and  VL(t)  is 

less  than  a  specified  tolerance  value. 

A  similar  voltage  update  scheme  can  be  used  to  improve  the 
multiple  reflection  technique.  Let  Fj(Vjj(t),V[j(t))  and 

FL<Vt>.  vu(t)}  be  the  reflected  voltage  waves  at  the  input  and 
output  terminals  respectively  in  the  original  multiple  reflection 


algorithm  (as  defined  In  section  3.2).  Using  the  voltage  update 
scheme 

Virk(maj)=Virk“,<mo>)+Pm(F|{Vi1k<t),VLik(t)|-Virk',(mco))  (3.15) 

VLrk(ma)=  VLrk' 1  (mo))+pm(FLtV1|k(t),VL1k(t))- VLrK~ 1  (ma>)) 
0.<Pm<!.,  m*0. 1,2,3, ... 

with  this  modification,  the  speed  of  convergence  can  be  increased 
significantly.  The  values  of  Pm  affect  the  speed  of  convergence.  The 

multiple  reflection  technique  is  the  limiting  case  of  this  algorithm 
with  Pm=1.  Typically,  the  optimum  values  of  Pm  is  .5  for  MESFET 
analysis. 


Chapter  4  Results 


The  validity  of  the  large  signal  model  is  confirmed  in 
two  ways.  First,  a  small  signal  equivalent  circuit  at  an  arbitrary  bias 
point  can  be  derived  from  the  large  signal  model  The  5-parameters 
of  this  small  signal  equivalent  circuit  are  calculated  and  compared 
with  the  measured  S-parameters(supplied  by  Hughes  Air  Craft)  of  the 
device  at  that  bias  point.  Second,  the  large  signal  model  is  simulated 
in  an  amplifier  confuguration.  The  power  input  versus  power  output 
curve  is  simulated  and  is  compared  with  measured  datavsupplied  by 
Hughes  Air  Craft).  Using  the  methods  in  chapter  2,  the  circuit 
elements  and  parameters  are  obtained  and  are  listed  in  table  4.1. 

4. 1  5-Parameters  comparison 

The  parameters  Gm  and  R0  in  the  small  signal  equivalent 

circuit(figure  2  1)  can  be  derived  from  the  large  signal  model  at  a 
bias  condition 


Gm=  *  'ch^g.Vd) 
d  Vg 


VGG-  VDD 


(4  1) 


R0=d  lch(Vq,Vd) 
d  Vd 


where  VGG  and  VDD  are  the  gate  and  drain  bias  voltages.respectiveiy. 

From  the  linear  element  values  in  table  41  and  Gm  ,  R0  as  derived, 

the  S-parameters  of  this  equivalent  circuit  can  be  calculated 
conveniently  using  Super-Compact.  The  S-parameters  thus  calculated 
for  the  PF-6000  MESFET  are  shown  in  figure  4.1  along  with  the 
measured  S-parameters  of  the  device  at  VGG— 2.03V,  VDD-7.V. 

4.2  Large  signal  simulation 

The  validity  of  the  modified  multiple  reflection  method 
is  first  verified  by  comparing  it  with  the  MESFET  simulation  results 
of  Materka  and  Kacprzak  [10].  The  circuit  parameters  of  Materkas 
amplifier  are  given  in  table  4.2.  The  model  is  simulated  in  an 
amplifier  configuration  as  shown  in  figure  4.2.  The  termination 
impedances  are  given  in  table  4.3.  The  results  are  shown  in  figure 
4.3.  By  means  of  the  modified  multiple  reflection  method,  the 
solution  to  this  particular  simulation  usually  coverges  within  5 
iterations.  When  the  original  multiple  reflection  method  was  used, 
the  soultions  converged  after  18  iterations. 

The  model  of  PF-6000  MESFET  is  also  simulated  in  the 
amplifier  configuration.  The  input  has  a  50  ohm  termination  and  the 
output  has  a  ll*8.5j  termination  impedance  at  15  GHZ.  Figure  4.4 
shows  the  simulated  and  measured  output  power  versus  input  power 
curves.  The  simulation  takes  an  average  of  7  iterations  to  converge 
with  the  modified  multiple  reflection  method.  When  the  original 


Table  4.1 


Model  parameters  for  Hughes'  PF-6000  MESFET 


Rs  .75  Q 

cds  W 

9o  279 

lsr  13mA 

Rd  80Q 

Csp  079pf 

q  2.56 

p  231 

Rg  1.14Q 

Cso  -613pf 

Vp0  -5.28v 

1  so 

Rj  .65  Q 

Cgo  -645pf 

y  -227 

ot  34 

Ls  .0 1 2nH 

Cdp  079pf 

Vs  2.23 

Ld  ,092nH 

*dss  496mA 

Vsat  1.85 

Lg  .  1 27nH 

c  4.52 

t  4.7ps 

Table  4.2 

Model  parameters  for  Materkas  MESFET 


Rs  4.5Q 

cds  ipf 

g0 

lsr  65nA 

Rd  4.5q 

^sp 

q  - 

p  1.28 

Rg  4.5Q 

^SO 

Vpo  -1.78V 

'so  145nA 

Rj  10.Q 

cgo 

g  -.11 

a  23 

Ls  InH 

Cdp  026pf 

Vg 

Ld  ,2nH 

1  dss  7^mA 

< 

CO 

cu 

r-f 

1 

1 

1 

1 

1 

Lg  ,2nH 

c  3.35 

T  5,ps 

Table  4.3 

Termination  impedances  for  Materka's  amplifier  simulation 

VGS=-75v  Vgs=-!.25v 


Pi(mw) 

hi 

-ss 

4l 

5. 

1 6.8+32.3  j 

67.9+67.6J 

24.9+ 

15.8J 

55.0+85.9J 

10. 

1 8.4+33.  Ij 

62.7+41. 9J 

38.8+ 

17  7j 

67.6+62.7j 

15. 

19.3+32.7j 

60.2+27.3J 

48.3+ 

18.7J 

68. 9+45. 6j 

20. 

19.8+32  5j 

55.8+19.8J 

55.8+ 

19  1  j 

64.9+35.2J 

25. 

20.0+3 1.7J 

53.6+ 1 6.  lj 

62.1  + 

19.4J 

61.1 +27.0  j 

30. 

20.2+3 1.4j 

49.3+13.0J 

67.5+ 

1 9.6  j 

57.9+22.7j 

35. 

21.4+31.5J 

47.2+11. 3  j 

72.2+ 

19.7J 

56.4+1 8  Oj 

40 

23.0+30.9J 

47.5+  10.5J 

76.1  + 

19.6J 

54  1  + 1 4.9j 

Table  4.3a  VD5=4.V 


VDS-30V  Vds-5.0V 


Pj(mw) 

hs. 

4l 

hs. 

4l 

5. 

18.6+32.8J 

68.1* 

44.3  j 

15.8+31  8J 

59.0+848J 

10. 

19.3+33.6J 

54.4+24.5J 

1 8.4+32  OJ 

68.4+62.7J 

15. 

19.9+33.3J 

47.6+ 

14.5j 

19  7+32  3j 

70.8+42.7j 

20. 

20.3+32.8J 

44.2+ 

9.6J 

20. 1  +32.4J 

640+32.3j 

25. 

20.3+31. 2  j 

42.3+ 

9  1 J 

20.6+31  8J 

62.0+26  6j 

30. 

20  9+31  3J 

41.4+ 

59J 

20.7+3 1.1  J 

60.0+20  7j 

35. 

22.3+30.8J 

41.0+ 

5.2J 

21.5+30.6J 

59.3+ 1 6.8  j 

40. 

23.9+40. 7J 

407+ 

4.9J 

22.7+30  2J 

59.1  +  15.8J 

Table  4.3b  vr. c=-.75V 


Po 

(mw) 


U  V 


—  Materka's  result 

—  simulated  result 


0  10  20  30  40 


Pi  (mw) 

Fig  4.3a  Output  power  vs.  input  power  of  Materka's  amplifier, 
f=9. 5  GHz,  Vds=4V. 


Vds=5V 


Pi  (mw) 

Fig.  4. 3b  Output  vs.  input  power  of  Materka's  amplifier, 
f=9 . 5  GHz,  Vgs=- .  75V 
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multiple  reflection  method  was  used,  the  solution  still  does  not 
converge  after  50  iterations,  even  for  a  small  input  power. 

When  simulating  the  circuit  of  liaterka  and  Kacprzak,  the 
termination  impedances  are  matched  to  the  FET.  Multiple  reflection 
method  works  fairly  well  when  termination  impedances  are  matched 
In  the  second  case,  the  input  termination  to  the  PF-6000  MESFET  is 
50  ohm  and  is  mismatched,  the  multiple  reflection  method  runs  into  a 
convergence  problem.  The  modified  multiple  reflection  method  shows 
good  covergence  speed  in  both  cases. 


Chapter  5  Conclusion  and  suggestion  for  further  research 


The  modeling  method  in  chapter  2  is  shown  to  be 
effective,  yet  relatively  simple  to  implement.  This  approach  is  very 
general;  it  can  be  applied  to  MESFETs  with  various  dimensions  and 
physical  parameters.  A  powerful  computer  aided  design  system  can 
be  developed  by  automating  the  circuit  element  and  parameter 
extracting  procedures  and  linking  then  together  with  Super-Compact 
and  the  curcuit  sumulator. 

The  modified  multiple  reflection  method  is  shown  to  be 
efficient.  The  accuracy  of  this  method  is  confirmed  by  comparing  its 
results  with  the  experimental  data  and  the  simulated  data  of  Materka 
and  Kacprza k.  The  validity  of  the  large  signal  model  is  confirmed  by 
comparing  the  simulated  and  measured  output  power  versus  input 
power  curves  of  the  PF-6000  MESFET.  On  the  average,  there  is  about 
3  dBm  difference  between  the  measured  data  and  sumulated  results 
This  difference  is  expected  to  be  contributed  by  3  factors: 

1.  Inadequacy  of  the  model. 

2.  Errors  involved  in  circuit  parameter  and  element  extracting 
procedures. 

3.  Errors  in  output  power  versus  input  power  measurements 
and  in  load  impedance  measurements. 

As  stated  in  chapter  2,  every  circuit  element  of  the  model  exhibits 
some  degree  of  nonlinearity  in  practice.  In  the  present  model,  the 


drain  to  sourch  capacitance  Cds  Is  considered  as  linear  element  To 

improve  the  accuracy,  C^g  probably  should  be  modeled  as  a  nonlinear 

element.  C^g  constitutes  the  capacitance  due  to  the  drain  and  source 

electrodes  and  the  capacitance  due  to  the  Gunn  domain  dipole.  The 
Gunn  domain  capacitance  is  certainly  voltage  dependent  in  nature. 


5.2  Suaaestion  for  further  research 


Several  areas  of  further  research  related  to  this  topic 
can  be  suggested.  In  the  theoretical  arena,  the  Gunn  domain 
phenomenon  of  the  GaAs  MESFET  needs  to  be  investigated.  This  Gunn 
domain  effect  is  expected  to  be  one  of  the  limitations  of  the  GaAs 
MESFET.  Shur[5]  and  Sing[6]  have  developed  an  analytical  model  for 
the  Gunn  domain.  However,  their  model  has  not  been  validated  under 
large-signal  condition. 

Using  the  large  signal  model  developed  for  the  GaAs 
MESFET,  other  microwave  devices  containing  this  MESFET  can  be 
analyzed.  The  modified  multiple  reflection  technique  can  be  applied 
directly  to  MESFET  multiplier  analysis.  For  MESFET  mixer  analysis,  a 
two  steps  large-signal  local  oscillator,  small-signal  RF  signal 
approach  taken.  The  multiple  reflection  mehtod  can  be  used  to 
analyze  the  large  signal  MESFET  curcuit  with  a  local  oscillator  input 
power.  A  subsequent  small-signal  analysis  is  needed  to  calculate  the 
RF-IF  signal  coversion  characteristic  of  the  mixer. 

For  MESFET  oscillator  simulation,  the  harmonic  balance 


method  is  the  more  appropriate  technique  to  use.  This  is  due  to  the 
fact  that  the  oscillating  frequency  and  the  termination  impedances 
can  be  taken  as  variables  to  be  optimized.  The  optimization  algorithm 
proposed  by  Walt[24]  for  diode  analysis  is  shown  to  be  very  fast  and 
efficient.  Thus,  this  algorithm  can  be  exploited  for  oscillator 
analysis 

Another  interesting  topic  is  the  modeling  and  simulation 
of  a  dual-gate  MESFET.  Dual-gate  MESFET  can  be  utilized  in  mixer 
design  in  which  the  local  oscillator  signal  and  the  RF  signal  are 
applied  separately  to  the  two  different  gates.  This  results  in  a 
simpler  design  of  a  mixer  circuit.  The  dual  gate  MESFET  can  also  be 
used  for  a  self-oscilating  mixer  in  which  one  of  the  gate  is 
terminated  by  variable  short  or  open  circuit  to  produce  local 
oscillator  power. 
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